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Abstract 

Means are used in several applications from electronic engeneering to 
information theory, however there is no general theorem on how to extend 
a given M(x,y) mean function to multiple variable forms. In this article 
we would like to present a theorem, which gives one possible solution for 
this problem, for every M(x, y) mean function, acting on positive numbers 
and operators. 


1 Introduction 

The use of mean functions falls back to the early ages of mathematics, but they 
are used widely nowdays aswell. For example they have great importance in 
statistics, but also in volume visualization, in imaging surgery. 

The study of electrical network connections implied the introduction of paral¬ 
lel sum of two positive semidefinite matrices in [2]. Formerly, Anderson defined 
a matrix operation, called shorted operation to a subspace, for each positive 
semidefinite matrix. Anderson and Trapp in [3], have extended the theory of 
parallel addition and shorted operation to bounded linear positive operators on 
a Hilbert space and demonstrated its importance in operator theory. They have 
studied fundamental properties of these operations. 

The axiomatic theory of means, for pair of positive operators, have been 
developed by Kubo and Ando in [7]. This theory has found a number of appli¬ 
cations in operator theory. 

Several steps have been taken in the theory of means, but a general idea has 
not been laid down yet, on how to extend, or define an n variable version of 
a given general M(x,y) mean function acting on numbers and operators. We 
would like to present a theory which gives a possible solution and a frame theory 
for further studies. 

Definition 1.1. A two variable function M : R + x R+ i— > R+ according to [11], 
is called a mean function if 
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(i) M(x , x) = x for every a; G R + . 

(ii) M(x,y ) = M(y, x) for every x,y G R + . 

(iii) If x < y, then x < M(x, y) < y. 

(iv) If x < x' and y < y', then M(x,y ) < M(x', y'). 

(v) M(x,y) is continuous. 

The geometric mean ^/xy, the arithmetic mean (x + y )/2 and the harmonic 
mean 2/(x _1 + y -1 ) are the most known examples but there are many other 
means aswell. 

The above definition of means can be easily extended to positive ordered 
operators and matrices, by replacing the numbers with them. 

2 Extending a mean to multiple variables 

The definition of M (x, y) given in the first section, can be extended to multiple 
variable functions. The following definition is one possibility of extension. 

Definition 2.1. An n variable function M n : (R + ) n i—> R+ may be called a (n 
variable) mean function if 

(i’) M n (x,... x) = x for every x G R + . 

(ii’) M n is independent from the ordering of xi, X 2 , ■ ■ ■ x n . 

(iii’) min(xi ,... x n ) < M n (x 1 ,... x n ) < max{x \,... x n ). 

(iv’) If 3(i) Xi < x\, then M n (x 1 ,... x„) < M n (x 1 ,... x\, ... x n ). 

(v’) If V(i) Xi <x\, then M n (xi,...x n ) < M n (xi,... x' n ). 

(vi’) M n (x 1 ,... x n ) is continuous. 

Our goal is to algorithmically define an n variable mean by using its less 
than n variable forms. Firstly we will use the n — 1 variable form of M, to 
define the n variable one, as an iteration’s limit. 

Definition 2.2. Let X = (xi,... x n ) G (R + ) n and M n _ 1 be an n — 1 variable 
mean function. Let us consider the n — 1 class variations of x n . There are 
( "jJ = n different variations. Let us define an iteration as 

x^ = Xi Vz G [1,... n] (2.1) 

xf +1 = \l„ ..t(V*(A' A )) (2.2) 

where V l n ^ 1 (X k ) = V^ n ^ 1 (x^,... x k ) is the zth n — 1 different variation of x n 
(n — 1 are chosen from n, which could be done in n different ways). 
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Theorem 2.1. The x k sequences - defined in definition 2.2 - are convergent 
and their limits are the same, which could be defined as the M n mean of the x n 
numbers. 


Proof. The iteration given in definition 2.2 has a contractive-like property by 
(iii’), which means that the sequence minX fc is monotonic increasing, rna xX k 
is monotonic decreasing. Hence the limits lim^oo min X k and lim^oo max 
exist. From condition (iv’), one can see that the series’ minimal and maximal 
elements are given by the following: 


min X k+l = M n ^(V^(X k )) (2.3) 

maxl t+1 = M n _i(V£-x (X k )) (2.4) 

where V™f n (X k ) are the smallest n — 1 numbers, V£-£(X k ) are the largest 
n—1 numbers from the series X k . By (2.3) and (2.4) minX fc+1 and maxX fc+1 
explicit dependence on min X k and max X k is given 

lim nxmX k+1 = lim M n _x(V^(X k )) (2.5) 

k —>-oo k—> oo 

lim maxl t+1 = lim M n _x(V^{X k )) (2.6) 


and by (vi’) 


one can write 

lim minX fc+1 

k—> oo 

lim maxl t+1 

k—> oo 


— M n _i 

— M n _i 


(lim VZ~!(X k )) 

k—>oo 


( lim V^{X k )) 

k—>oo 


which yields 


lim minX fc+1 = / = M n _i(/,...) 

k—> oo 

lim maxl i+1 = J = M„_i(J,...) 


(2.7) 

( 2 . 8 ) 


(2.9) 

( 2 . 10 ) 


and this is only true when I = J. 


□ 


The above theorem and its proof yields the following remarks. 

Remark 2.1. The iteratively defined mean function M n , is invariant to the initial 
ordering of the V) n_1 variations. 

Remark 2.2. The iteration (2.2) leaves the mean of the starting n numbers 
invariant through the sequence. 

It is easy to verify the next two theorems, which have stressed importance 
in inequalities of means and operator means in our given context. 


Corollary 2.2. Two different, iteratively defined mean function M n i and M n 2 , 
is in the same relation as their two variable forms (Mi x < ^ 2,2 implies M n \ < 

Mn t i). 
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Corollary 2.3. Theorem 2.1 and the proof also works for ordered positive op¬ 
erators, acting on a H Hilbej't space, and for k x k matrices. 

We will show with some examples, that theorem 2.1 gives a sufficient defi¬ 
nition of the n variable mean. 


Corollary 2.4. Theorem 2.1 applied on the n— 1 variable arithmetic, geometric 
and harmonic mean, gives the corresponding n variable mean. 


Proof. According to the given serie’s convergence in theorem 2.1 it is enough to 
prove that the minimum’s or maximum’s limit is the n variable mean. We will 
prove it only for the arithmetic mean. Let us consider the numbers x'( < ■ ■ ■ < 
£ K + and the A„_i arithmetic mean. By theorem 2.1, the sequence x\ can 
explicitly be written and proven by induction with (2.2), (2.3) and (2.4): 


x 


k 

l 


x 


k 

n 


E?=i- 

(n-l) fe -l 0 | 0 


Liu 

(n—l) fc 

(T ,_l)fc + l ^ 0 

y-v n 

M i=1 

(n—l) fc 

(n-l)fc-l 0 0 

n i '' tj n 

£?=i- 

(n- l) k 

(n — l)' t +l x 0 ^,0 


(n — 1) h 


if k is even, 
if k is odd, 

if k is even, 
if k is odd. 


( 2 . 11 ) 


For the geometric mean the proof can be extended using the logarithmic 
function and its inverse for the limit: 


log( n f/x i • x 2 ■ ■ -x n -i) = 


T,?=i log^ 


n — 1 

The proof for the harmonic mean can be given by inverses: 

-l 


n — 1 


-i , -l 


Xr, ... 


x 


1 Xcy 1 


. X 


n —1 


n — 1 


( 2 . 12 ) 


(2.13) 

□ 


Our main idea of extending means to multiple variables is based on theorem 
2.1, which is theoretically enough but in practice is very insufficient. For example 
if we would like to compute n numbers or matrices M n mean , we should use the 
two variable main definition M 2 and extend the other ones from one to another. 
In the next section we will prove that M n can be directly extended from the 
corresponding M 2 - 


3 Extending M n directly from M 2 

Let us define the following iteration: 
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Definition 3.1. Let X = (x\ < ••• < x°) £ (R + ) n and M = M 2 be a two 
variable mean function, 

2 ) if i = 1, 

M{ x n-n x n) if i = n, (3.1) 

M[x\_ i,4 + i) else. 

Theorem 3.1. T/ie iteration given in definition 3.1 for all n is convergent and 
V(i) lim a;* = M n (x°,... x° n ) 

k—>00 


where M n is defined by theorem 2.1. 

Proof. Firstly we begin with proving the convergence. It is clear that the min x k 
is always the first element (i = 1) and the maxxf is always the last element 
(i = n) of the series in definition 3.1. Hence (iii) and definition 3.1, mina:^ is 
increasing and max x\ is decreasing. This yields: 


and by (v): 


which give 


lim mina;^ +1 = lim M(x k , x k ) 

k—>oo k—>00 

lim maxij +1 = lim M( x k _ 1: x k ) 

k —too k —>00 


lim mmXi +1 = M( lim x k , lim x k ) 

k—>00 k—>00 k—>00 

lim nraxa:^ +1 = M( lim x k _ x , lim x k ) 

k —too k—>00 k —too 


lim mina;* +1 = I = M(I, lim x k ), 

k—>oo k—>oo 

lim maxa;^ +1 = J = M{ lim x^_ 1: J ). 

k—>oo k—>00 


(3.2) 

(3.3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 


Considering the characteristics of the definition 3.1, this can only be true when 
I = J. 

Secondly we will prove the limit. For n = 3 the theorem is clear, because 
the two iterations, defined in theorems 2.1 and in 3.1, are the same. Our next 
step is to prove for n + 1, if it is true for n. 

Let us consider the definition of M n+ 1 in theorem 2.1. Comparing the 
rnina:^ (which is the first element i = 1) in theorem 3.1, and min X k (which 
equals M n (Vf l l lin (X k ~ 1 )) by (2.3)), we can see that mina;-' < minX fe , be¬ 
cause of the inductional condition and the definition of M n as a limit in the¬ 
orem 2.1. The same can be applied for maxa:^ and ma xX k which yields 
maxa:f > max X k . Hence mina:f and maxa:^ are minoring and majoring, 
for every k, mm.X k and ma xX k , but lim^^ min x k = lirn^oo max x k , so 
limfc^oo min x k = M n+ i(a^,... a;° +1 ). □ 
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Furthermore there is special property in the iteration in definition 3.1. 

Definition 3.2. Let x'f € R + i € [1,... n] and G be a graph, with n verteces 
and edges given as that, there is one cycle in G, which contains all verteces and 
edges (so it is at the same time a Hamiltonian- and an Euler-cycle). This implies 
that in G, every vertex has two edges and all of them are bound together. Let 
us consider an optional one to one correspondence between x® numbers and G-s 
verteces. Taking every edge in G as an M(xj, xf) (where M is a mean function 
and Xj, x ° are assigned to the two ending points of the edge as previously given), 
we can define an iteration with an optional n mappings, 


a* +1 =M(x*,xt) 
i,j,lG [l,...n] j^l. 


(3.8) 


Theorem 3.2. Every different iteration given in definition 3.2, converge to the 
limit M n mean function, and the iteration - independently from the mapping 
x^ +1 = M(xj, xf) - converge on a higher or equally rate as the iteration given 
in definition 3.1. 

Proof. For n = 3, it is easy to see that the theorem is true, because the iterations 
given in definitions 3.1 and 3.2, are the same. 

Assume that the theorem is true for n variable. Let us expand from an n 
variable iteration defined in 3.2 with an optional mapping, to n + 1 variable. 
This can be done as replacing one edge with two edges and one vertex (mapped 
to a new number). Let us do this expansion in the following way. Take the 
first smallest n numbers from n + 1 and set up on them the iteration given 
in definition 3.1. From the inductional condition this iteration will have the 
slowest convergence rate, which means that its minimal and maximal elements 
will minor and major every other iteration in definition 3.2. Let us replace 
the edge which gives the maximal element of the iteration given in definition 
3.1, with the two new edges and the remaining number (which is the greatest 
number out of the n + 1) as a vertex. This two edges with the corresponding 
two M 2 (x,y) will give the new iterations - given for n + 1 numbers - greatest 
two elements. This replacement cannot be done better in any other optional 
mapped iteration given in definition 3.2 aswell. But considering the inductional 
condition this yields that any n + 1 variable iteration given in definition 3.2, 
cannot minor and major, with its maximal and minimal elements, the iteration 
given in definition 3.1, hence theorem 3.2 is proven. □ 

Corollary 3.3. The above theorems also work for ordered operators acting on 
a H Hilbert space, and k x k matrices aswell. 

We will have to consider further examinations to define the above itera- 
tional definitions for inorderable matrices. In the next section we will study this 
problem. 


6 



4 Extending theorems 3.1 and 3.2 to unordered 
matrices and operators 

The problem is with positive matrices and operators which satisfy ||A|| = ||.B|| 
and A ^ B. For the above matrices and operators, the function M(A,B)’s 
and its arguments’ relation is not explained and highly depend on the main 
characteristics of M(A,B), so the given iterations in the above theorems must 
be specified. 

Theorem 4.1. For any X±,... X n positive operators or matrices the iteration 
given in definition 3.1 is convergent, as defined in theorem 3.1. 

Proof. If ||Xi|| = ||X 2 || = ••• = ||X„|| does not hold, than the iteration in 
definition 3.1 converges for all Xj, because after n steps, the iteration will surely 
alter all of the ||X^||-s (from one to another), so the iteration will converge. 

If ||Xl|| = ||X 2 II = • • • = 11 X n 11 does hold, we will have to define (according 
to [11]) the following construction. Let a t and a' t be monotone sequences as, 

at,a t € R + 

V(t) a t > 1 and a' t < 1 (4.1) 

lim t ^oo a t = 1 and lim t ^= 1. 

Let X[ = a t X u X[ = X, i e [2,... n] and X'{ = a' t X 1 , Xf = X t i e [2,... n\. 
Let us set the iteration given in definition 3.1 up on X it X- and X". According 
to the first part of the proof, the (X') fc and (X") fc series are convergent for any 
t, as given in theorem 3.1. Considering the definition of sequences a t and a' t 
in (4.1), it is easy to verify by condition (iii), that for any t, the (X") fc series 
are minoring and {X'f) k series are majoring the series X k for any i. Taking 
the limit k — > 00 , we get M n (X"(t),... X"(t)) and M n (X[(t),... X' n (t)). Hence 
M n (X"(t),... X"(f)) and M n (X[(t),... X^(t)) are Cauchy sequences in index t 
and condition (vi’), they are convergent and 

lim t ->ooM n (Xi(t),... X"(t)) = limt^ooM n (X[(t),... X' n {t)). (4.2) 

But (Xf) k (t) and (X[) k (t) are minoring and majoring every X k for any i and 
t, so the limit limk->ooX k exist and by (4.2), 

lim t ^ooM n (X"(t),... X"(t)) = Um.k^ooX k = lim t ^ooM n (X[{t),... X' n (t)) 

(4.3) 

and theorem 4.1 is proven. □ 

Corollary 4.2. Using the above proof, theorems 2.1 and 3.2 work for the un¬ 
ordered Xi-s. 

5 Consequences 

By the theorems given in our examinations generalize the extension of the two 
variable mean functions and gives a frame theory, which may be used in the 
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future studies related to the extension of means to multiple variables. An im¬ 
portant outcome is, that these theorems are applying for operators and matrices 
and guarantee the existence of one possible extension. 

It is known, that in several situations, there are more than one possible 
generalization of a mean. One example is the logarithmic mean, 


L(x,y) 


x-y 

log x — log y 


(5.1) 


which has several extended forms according to [5], [8], [9], but our theorems may 
leave only one form valid. However with some means , it appears to be quite 
difficult to give the iterations limit in a closed form. 
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